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The dynamics of an incompressible, dissipationless Hall magnetohydrodynamic medium are in¬ 
vestigated from Lagrangian mechanical viewpoint. The hybrid and magnetic helicities are shown to 
emerge, respectively, from the application of the particle relabeling symmetry for ion and electron 
flows to Noether’s first theorem, while the constant of motion associated with the theorem is gen¬ 
erally given by their arbitrary linear combination. Furthermore, integral path variation associated 
with the invariant action is expressed by the operation of an integro-differential operator on the ref¬ 
erence path. The eigenfunctions of this operator are double Beltrami flows, i.e. force-free stationary 
solutions to the equation of motion and provide a family of orthogonal function bases that yields 
the spectral representation of the equation of motion with a remarkably simple form. Among the 
double Beltrami flows, considering the influence of a uniform background magnetic field and the 
Hall term effect vanishing limit, the generalized Elsasser variables are found to be the most suitable 
for avoiding problems with singularities in the standard magnetohydrodynamic limit. 

PACS numbers: 52.30.-q,45.20.-d,52.35.Mw,47.10.-g 


I. INTRODUCTION 

In the present study we investigate dynamical system 
features of a dissipationless incompressible Hall magneto¬ 
hydrodynamic (HMHD) medium and propose the notion 
of helicity-based, particle-relabeling operator, which is lo¬ 
cated at the junction of two seemingly separated topics: 
particle relabeling symmetry and force-free, stationary 
state solution. Consideration of the invariant action as¬ 
sociated with the particle-relabeling symmetry naturally 
leads to the operator, and its eigenvalue problem and as¬ 
sociated normal-mode expansion of basic formulas and 
equations are examined as its application. 

The HMHD is well-known as a simple, one-fluid plasma 
model that contains two-fluid effects and that has been 
intensively investigated both numerically and mathemat¬ 
ically. The basic idea of the HMHD approximation is for¬ 
mulated by replacing the magnetohydrodynamic (MHD) 
approximation [lj, by which a vanishing Lorentz force is 
assumed for the entire plasma {E + V x B = 0), with 
an assumption that the Lorentz force only vanishes for 
the electron component of the plasma (E + V e x B = 0), 
where V and V e are the averaged plasma velocity and 
its electron component, respectively !2j. The formula¬ 
tion is completed by approximating the entire plasma 
velocity by its ion component (V ss Vf), and then eval¬ 
uating the current density by V) — V e = J/en e with an 
approximated Ampere’s law, V x B = RoJ- Thus, the 
evolution equations for a dissipationless, incompressible 
HMHD plasma are given by the incompressibility con¬ 
dition, solenoidal condition, momentum equation, and 
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induction equation as follows: 

V • u = V • b = 0, 

dtu = u x (V x u) + j x b — VP, (1) 

dtb = V x ((m — aj) x b ), 

where u , b , j, P, and a are the appropriately nondimen- 
sionalized variables corresponding to ion velocity, mag¬ 
netic field, current density (j = V x 6), generalized pres¬ 
sure, and Hall term strength parameter, respectively. It 
is easy to see that, in the limit a —>• 0, the system reduces 
to the standard MHD system. 

The system HD is known to have three constants of 
motion, i.e., the total energy, E, the magnetic helicity, 
Hm, and the hybrid helicity, Hjj, which are given by 
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respectively Q, where a is the vector potential of b (b = 
V x a). Obviously, in the MHD limit, a —s- 0, the hybrid 
helicity degenerates into the magnetic helicity. However, 
it was very interesting that the spectral representation 
of Q]) by the generalized Elsasser variables was naturally 
proved to yield four constants of motion due to the skew- 
symmetry of the quadratic terms coefficients Q. The 
fourth constant was the modified cross helicity, given by 
He ■= Hh — Hm, and it converges to the cross helicity 
in the MHD limit. Despite that the helicity conservation 
has been known to emerge from the particle relabeling 
symmetry for the MHD case ||, its HMHD counterpart 
still remains unresolved. 
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We focus here on the Lagrangian mechanical aspects of 
the HMHD system together with a differential topological 
framework. 

Since Holm established the Hamiltonian mechanical 
description of the HMHD system Q, analytical mechan¬ 
ical approaches to HMHD physics have been mainly car¬ 
ried out within the Hamiltonian mechanics framework 
[ 3 , 1 - Recently, a Lagrangian mechanical approach was 
employed by Keramidas Charidakos et al. [§]. Their 
Lagrangian was obtained by naturally extending an 71- 
particle system Lagrangian to a two-fluid plasma model, 
and the HMHD momentum equation and Ohm’s law were 
derived. 


On the other hand, our Lagrangian mechanical ap¬ 
proach is rather close to Arnold’s differential-geometrical 
method. 

Ever since Arnold [Toj ] reviewed his studies of dynami¬ 
cal systems on Lie groups and related hydrodynamic top¬ 
ics in a unified form, many fluid dynamical systems have 
been recognized to exist on appropriate Lie groups Q- 
The key mathematical objects of Arnold’s method are 
twofold. One is an appropriate Riemannian metric that 
is introduced on the relevant Lie group as a Lagrangian 
of the action. The other are so-called “Lin’s constraints” 
that provide the variation of an integration path [ 12j . 

In the field of plasma physics, Arnold’s method was 
found to be applicable to the dynamics of a dissipa¬ 
tionless, incompressible MHD medium if the Lie alge¬ 
braic structure was appropriately defined on the func¬ 
tion space of the pair of the velocity and magnetic fields 


[131114| . This extension to the pair was called “magnetic 
extension” and is now recognized as a special case of the 
semidirect product of a Lie group and a certain vector 
space (Sect. 10.B of Ref. 11|). 

Since the induction equation is not “passive” due to 
the Hall term (i.e., the magnetic field can evolve au¬ 
tonomously) , the HMHD system does not obey the mag¬ 
netic extension scheme. This discrepancy of magnetic 
extension can be overcome by replacing the group ac¬ 
tion on the vector space with the group homomorphism, 
which was based on Vizman’s extended formulation fl5| . 
The configuration space was given by a semidirect prod¬ 
uct of two volume-preserving diffeomorphisms and the 
Lagrangian was given by a Riemannian metric that 
physically implied the total plasma energy Q. In the 
present study, we will report another HMHD formula¬ 
tion, wherein the ion and electron velocities are taken 
as basic variables and the configuration space is given 
by a direct product of two volume-preserving diffeomor¬ 
phisms, and we will discuss the conservation of hclicities 
as a consequence of the particle relabeling symmetry of 
each fluid. For the MHD case, the particle relabeling 
symmetry and its relation to helicity conservation is dis¬ 
cussed by Padliye and Morrison Q. 

Note that, despite the simple appearance of the ba¬ 
sic equations ©, analytical mechanical approaches to 
HMHD systems raise a small parameter problem when 
their relation to the standard MHD limit is considered. 


For example, in the Hamiltonian mechanics approach, 
one of the natural choices of vector variables is the pair 
of the total ion momentum density, M = pv + R~ 1 apA , 
and the magnetic vector potential, A, where p , v are the 
density and velocity of ion component and R/a = a in 
our notation ||. In the limit a —* 0, these two variables 
come close to each other, M R _1 apA, and manip¬ 
ulation of the small difference v = M/p — R~ 1 aA is 
needed to capture the ion flow. Recently, Yoshida and 
Hameiri proposed a method to treat the MHD limit by 
renormalizing the Lagrangian described by some appro¬ 
priate Clebsch variables [16j . In the present study, we will 
seek another way of avoiding the singularity problem by 
choosing an appropriate expansion function set. 

In the context of the analysis of fully-developed turbu¬ 
lence, it was recently shown by direct numerical simula¬ 
tion (DNS) that the Hall term effect alters the formation 
tendency of coherent structures m. Formation of tubu¬ 
lar structures of currents and enstrophy densities at small 
scales are observed for the HMHD case, while sheet-like 
structures are often observed for the standard MHD sys¬ 
tem. In addition, it is interesting that although both the 
Lorentz force term of the ion velocity evolution equa¬ 
tion and the Hall term of the magnetic field evolution 
equation contain the function j x b, their contributions 
to the energy transfer of the kinetic and magnetic ener¬ 
gies were found to be quite different [l8j]. This suggests 
that, for the analysis of basic dynamical features, it is 
not sufficient to focus upon the features of a magnetic 
field alone, but that their coupling with the velocity field 
must also be considered. Thus, an appropriate coupled 
base function system is required for the DNS or some 
other practical analysis. 

In relation to the coupling of the magnetic and ion 
velocity fields, there exist two significant functional cat¬ 
egories to describe the equilibrium states, dynamics, and 
the stability of the HMHD system: the double Bel¬ 
trami flow (DBF); and the generalized Elsasser variable 
(GEV). 

The notion of DBF was introduced by Mahajan and 
Yoshida in order to extend the concept of a Taylor state 
to two-fluid plasma models [19]. To derive the stable 
equilibrium state, the DBF was applied to the variational 
calculation to minimize a dissipation function j2Cj. The 
DBF was also applied as “dynamically accessible varia¬ 
tion,” which conserves the Casimir invariants, to analyze 
the nonlinear stability of the equilibrium state [§]. In 
the present study, the Casimir-preserving nature of the 
DBFs will be considered based on its Lagrangian me¬ 
chanical counterpart, i.e., Noether’s first theorem, and 
the eigenfunctions of the DBF-generating operator will 
be shown to provide a remarkably simple expression for 
the evolution equation. 

On the other hand, the GEV was introduced by Galtier 
to formulate the HMHD dynamics in the wave/weak tur¬ 
bulence closure analysis framework [21]. Though GEVs 
were developed to describe the linear waves that are ex¬ 
cited when a uniform background magnetic field exists, 
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they can also be used as a set of orthogonal base func¬ 
tions even when the ambient field is absent. In the pre¬ 
vious study, it was shown that the GEV expansion of the 
HMHD equation naturally yields four conservation laws 
due to the symmetric properties of the quadratic term 
coefficients, and it was conjectured that this observation 
might reflect some symmetries intrinsic to the system Q. 
Recently, we also have applied the GEV decomposition 
to the DNS data and confirmed the mirror symmetry 
breaking at small scales [22:]. In the present study, we 
will review the GEV as a specific example of DBF and 
discuss its advantages over other DBFs in relation to the 
MHD limit. 


This paper is organized as follows: the basics of the La- 
grangian mechanical formulation are given in section 2; 
variational calculations are carried out in section 3, where 
we derive the equation of motion from Hamilton’s princi¬ 
ple and the conservation of the helicities from Noether’s 
first theorem; in section 4, the derivation process is re¬ 
formulated using the differential topological terminology, 
and the topological foundations of hclicity conservation 
are discussed. The DBFs are used as the base functions of 
the HMHD system and the topological basic quantities, 
i.e., the Riemannian metric and the structure constant 
of the Lie group, are given in the section 5; the influence 
of a uniform background magnetic field and the standard 
MHD limit of HMHD system are discussed in section 6; 
the section 7 is devoted to discussing the implications of 
our findings. 


II. FORMULATION 


In a Lagrangian mechanical description of hydrody¬ 
namics, the basic variable for describing the fluid motion 
is known to be given by an n-tuple of functions that maps 
the fluid particles from one time to another; we call this 
variable the “particle trajectory map” (PTM) hereafter. 

In the present study, we choose as the basic variables, a 
pair of PTMs, say (X(t), Y(t)), which describe the posi¬ 
tions of the ion and electron fluid particles, respectively; 
X{a,t), Y(a,t) express the positions, which are initially 
(at t = 0) located at a € M.\ 23] In differential topologi¬ 
cal terminology, we consider here the dynamical system 
on a direct product of two volume-preserving diffeomor- 
phisms, say G := 5'Diff(M) x ST)ifF(M), hereafter. Note 
that, the choice of basic variables is not unique for the 
HMHD system; in our previous study, we used the pair 
of PTMs of the ion velocity and the current density to 
constitute a semidirect product of diffeomorphisms [4j. 

The PTMs are related to the ion and electron velocity 
fields in the Eulerian specification, say Vi (t) = V k ( t ) , 


Ve(t) = V k (t)° e £ S (M), by 
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hereafter, 3fs(M) denotes the function space of the 
divergence-free, tangent vector fields on M. Mathemat¬ 
ically, the RHS’s of these equations express the right 
translation of the vector field, V. (resp. V e ), by the group 
operation X(t) (resp. Y(t)). As was discussed in f24| . 
since the arguments of component function and basis do 
not agree with each other, the LHS’s of © are not proper 
differential topological objects; the Lagrangian velocities 
in the RHS’s form are appropriate for the calculus on 
manifolds. 

In Lagrangian mechanics on Lie groups, there exist 
two key mathematical structures: the Lie bracket and 
the Riemannian metric. The Lie bracket is necessary to 
determine the higher-order terms of the Taylor expan¬ 
sion of a composite function of PTMs. The Riemannian 
metric is the inner product of two tangential vectors of 
PTMs and defines the Lagrangian of the system. 

Since the group operation of G is defined by the 
compositions of function triplets (Xi. Y\) o (X 2 , Yi) = 
(Xl(X 2 ), Yi(l 2 )), the Lie bracket of the associated Lie 
algebra is given by 



(v x (V a x V i2 ),V x (V e i x V e2 )), 


( 6 ) 


where V k = (V lk ,V ek ) £ g = T e G = X S (M) x £ S (M). 
Since M is three-dimensional and the vector fields consid¬ 
ered here are divergence free, the Lie bracket on £s(M) 
is given by [a, b] = V x (a x b). 

In the present study, the Riemannian metric at 
(X, Y) € G is defined by the combination of the integrals 
described by the Lagrangian and the Eulerian specifica¬ 
tions as follows: 
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where d 3 X(a,t) and (Vx) -1 are the advected volume 
element at the time t (which is initially located at 
a) and the inverse of the curl operator, respectively. 
For practical calculations, the first term is replaced by 
L m d 3 f (Vj! • Va)-, because the modulus of the vol¬ 
ume element is conserved due to the incompressibility: 
|d 3 A(a, t)\ = |d 3 a| for all t. Mathematically, this re¬ 
placement implies the right invariance of the Riemannian 
metric. 

Since the difference V — V e gives a current den¬ 
sity aj, the generated magnetic field is given by b = 
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a _ 1 (Vx) _1 (Vi — V e \). Thus, the Riemannian metric 
expresses the sum of the kinetic energy of the ion flow 
(with density p = 1 ) and the magnetic field energy gen¬ 
erated by the plasma current, while the kinetic energy of 
the electron flow is assumed to be negligible. 

In the present formulation, the Riemannian metric co¬ 
ercively combines the two different vector spaces by a 
subtraction operation, although the implication of the 
operation is quite natural from a physical viewpoint. In 
our previous study, the coupling of two spaces is estab¬ 
lished by the group action of a semidirect product of two 
diffeomorphism groups, while the Riemannian metric is 
defined by separately defined integrals. 

A remark on the Lie algebraic structure should be 
made here. Substituting 

Vi = u, V e =u- aj (8) 


into d6j) and CD, we obtain the inner product of a In¬ 
variable and their Lie bracket as follows: 


[V 2 ,V 3 ]\ =[ d 3 4«r- [Vx (t * 2 xt* 3 )] 

' ( e > e ) JseM l 

+bi ■ (u'2 x j 3 + j‘2 x u 3 - aj 2 X j 3 ) | , (9) 


where bk satisfies bk = (Vx) -1 .?*., V ■ bk = 0. The same 
integral can be obtained from the other Riemannian met¬ 
ric and the commutator, i.e., from Eqs.(4) and (7) (or 
Eqs.(12) and (13)) of our previous study Q. This implies 
that these two formulations provide the same structure 
constant of the Lie algebra if appropriate base functions 
such as the GEVs are applied, and thus, these two sys¬ 
tems are equivalent although their group structures are 
quite different from each other. In other words, these 
two formulations constitute a kind of “canonical transfor¬ 
mation” between the configuration spaces with different 
group structures. 


III. VARIATIONAL CALCULATION: 
DERIVATION OF THE EQUATION OF MOTION 
AND HELICITY CONSERVATION 

Action along a path 7 (t) = (X(t),Y(t)) G G (t G [0,1]) 
is given by S = f 0 Ldt , where L is the Lagrangian defined 

by 

L = L{ 1 {t),i(t)):=\(V\V) iW (10) 

where j(t + r) r* (exp(rV(t)),exp(rV e (f))) o 7 (f), t G 
[0,1], t is a small parameter, and exp is the exponential 
map on £s(M). Let 7 (t;6) be a perturbed path, where 
7 M) =7 (t), 

7 (t;S) « (exp(^(t)),exp(<5r/(t))) o 7 (t;0), 

6 is a small parameter, and (£, rj) G 0 = Xs(M) xXs(M) 
are the displacement fields. Noticing that the perturba¬ 
tion part of the velocity, say (Vi,V e ), obeys Lin con¬ 
straints 


(see Appendix [A] for the derivation), the first variation 
of the action is given by 


SS = dt d J a?<J Vi • V 
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where A is the vector potential of the magnetic field with 
Coulomb gauge divided by a: 

A:=a~ 2 (V x)~\v~ V e ) = (13) 

In the present study, we assume that the boundary in¬ 
tegrals always vanish. The expression in the second line 
yields the following two results: first, the conjugate mo¬ 
menta of V, and V e are given by 


Mi := 

M e := 


SL 

6V 

SL 

6V e 


— Vi -}-A — u 


a 

a 
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respectively; second, the variations due to £ and r/ that 
satisfy 

^+Vx(^x V) =0, d t r 7 + V x (77 x V e ) = 0 (15) 

retain the value of action. In terms of Lin constraints 
CD, this reads as (V.,V e ) = (0,0), i.e., the velocity 
fields along the perturbed paths are the same as those 
of the reference path. This symmetry for the invariant 
action is well-known as the particle relabeling symmetry 
[IIH; we give a brief review in Appendix |B] 


By integration by parts of (fl2l) with respect to t and x 
and changing the order of scalar triple products of vector 
fields, we obtain 


SS = 


f dt f d 3 4 (v + A) ■ d t i + 

Jo Jx£M l 

4 • [Vi X (V x (v + A)) 

A • dtTj - 77 • [Ve X (Vx A)]J, 

f d 3 x ((V + A) • £ - A ■ rj) t=1 

J xGM 

-f d 3 x((V i + A)-Z-A.r 1 ) t=0 

JxGM 


Jx£M 

!>. 
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xU- — d t (Vi + A) + Vi x (V x (V + A)) 
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d t A - V e x (V x A) 


(V, Ve) = dt(£, r?) + [(e rj), (V, V e )] (11) 


(17) 
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Hamilton’s principle, i.e., SS = 0 for an arbitrary pertur¬ 
bation (£, r)) with fixed path end conditions, leads to the 
following Euler-Lagrange equations: 


Noether’s first theorem: 


SS = 


L 


d 3 a 


SeM 


(V, + A) • (V x (V, + A)) 


d t (V + A) = V x (V x (V+A)) - VP,, (18) 
dtA = Ve X (V X A) - V P e , (19) 


/ d 3 z 

IxeM 
-2 TT fi\ _ —2 ] 


^[(VJ + A)- (Vx (V* + A)) 

= ar 2 Hff(l) - (0) = 0. 


t =0 

( 22 ) 


where Pi and P e are the generalized pressures for each 
fluid. Substituting @ and m, and carrying out some 
calculations, we obtain the evolution equation ©■ Note 
that, in the limit a —> 0, we obtain the standard MHD 
equations, although the variable A = a/a diverges at 
0(a~ 1 ). 


Next, we consider specific perturbations (£, rj) that 
leave the value of action unchanged. The combinations 
of $ = Vi or V x (Vi + A) and rj = V e or V x A are the 
candidates for the invariant action, because they can¬ 
cel the cubic terms of (fT()l) . It is easy to see that the 
conservation of total energy © is obtained by setting 
{ZiOiWit)) = (Vi,V e ), which implies that the variation 
is taken in the direction of the path, i.e., the variation is 
associated with the time translation. 

By setting (£(t),r](t)) = (0, V x A(t)), the first varia¬ 
tion of action (usd becomes 

SS = f dt [ d 3 x(~A-(Vxd t A)), 

= - [ dt [ d 3 f((Vxi)-9 f A (20) 
Jo JxeM ' ' 

where the second line is obtained by the integration by 
parts with respect to x. If the path 7 (t) satisfies the 
Euler-Lagrange equation (IT9l) . substitution of the RHS 
of m into (HU), yields the vanishing first variation; i.e., 
SS = 0. Noticing that the identity 
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holds, the integration by parts with respect to t without 
fixed path end conditions results in the conservation of 
magnetic helicity ([3]) by Noether’s first theorem: 


SS = 


dt 


[ d -4 

IxeM 


|-A.(Vxi 


d 3 


SeM 


IxeM 


|A'(VxA 
d 3 f(A- (V x A)) 


t= 1 


= a~ 2 H M (l)~a- 2 H M (0)=0. 


( 21 ) 


Similarly, by setting (£(t),rj(t)) = (V x (V(t) + A(t)), 0), 
we obtain the conservation of hybrid helicity 0 by 


In summary, if the path j(t) locally satisfies the 
Euler-Lagrange equations m and m, the combina¬ 
tions of perturbations (£,rj) = CE(Vi,V e ) and (C)V x 
(V) + A),— C e V x A) retain the value of the action, where 
Ce, Ci and C e are arbitrary constants. The derivation 
process clearly shows that the magnetic (resp. hybrid) 
helicity is obtained by varying the integral path only on 
the !4-(resp. V)~)side of the configuration space; i.e., 
the magnetic and hybrid helicities are obtained by the 
relabeling of V e and Vi. 


IV. DIFFERENTIAL TOPOLOGICAL 
DESCRIPTION 


In order to obtain mathematical insight to these con¬ 
servation laws, we revisit the discussion above using the 
differential topological expressions. 

As a starting point, we notice that the Lie bracket is 
also expressed by the Lie derivative of a vector field; i.e., 


DP Dtf d 

= = (23) 


for r) € £(M). Using this relation, the first variation 
m can be rewritten as 


SS = J dt{ (m\ (dt + L Vi )£) + (M e \(dt + L Ve )rj )}, 


(24) 


where the underline denotes a differential 1 -form and the 
parenthesis is the inner product between a differential 1 - 
form and a vector field; ( * | * ) : Cli(M) x X(M) —X R, 
where (A\B) = f SeM AiB l d 3 x, A = Aidx 1 £ fli(M), 
and B = B l -J=P £ X(M). Integration by parts yields 

SS = 

{(M^) + (Me\rj)} t=i {(Ml^) + (Me\d)} t=o 

dt{((dt + L Vi )Mi ||) + ({dt + L Ve )Me |r?)}, (25) 

which corresponds to m ■ Under the fixed path end 
conditions, £ = ij = 0 for t = 0 and 1 , we obtain the 
Euler-Lagrange equation [26] 

(d t +L Vt )M 1 = -dPfi (d t + L Ve )Me = ~dPfi (26) 

where P[ and P' e are introduced to satisfy the divergence- 
free condition. These are the differential topological ex¬ 
pressions of m and (HU). Since the exterior differentia¬ 
tion, d , is commutative with the Lie derivative of differ¬ 
ential form (24|, |27j , using dd = 0, we obtain the exterior 
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derivative of (1261) as follows: 

(< d t + L Vz )dMi=0 , (d t + L v JdMe = 0. (27) 

Since we consider three dimensional space and 
divergence-free vector fields and differential forms here, 
there is a natural correspondence between the differential 
2-form and the vector field [Hj]. Here, we introduce the 
mapping [*] : f h(M) —> £(M), defined by 

(AM) :=/^A|, (28) 

where A G Di(M), and £ G II 2 (M). Using the relations 

= [rf Mr l, which are guaranteed by the 
divergence-free condition, we obtain 

(d t + L Vi )[dM j] = 0, (dt + L Vr )\dM e ] = 0. (29) 

These equations obey the invariant action conditions 
(USD- Thus, by substituting £ = aC'Jri Mi l = CflaV x 
m + 6 ) and 77 = aC e [rfM e ] = C e b into ESl) . we obtain 
the general helicity conservation law H(l) — 17(0) = 0 as 
a consequence of invariant action, where the constant H, 
which we call the mixed helicity hereafter, is given by 

H = aCi (Mi | \dMj}) + aC e (Me\ \dM e }) , (30) 

and Ci and C e are arbitrary constants. Using (l28l) . H 
is also written as the integral of the wedge product of 
differential forms: 

H = a f (CjMj A dMj + C e M e A dM e ). (31) 



FIG. 1. Relation among generalized velocity, V, helicity- 
based, particle-relabeling operator, W, and act ion-preserving 
path variation, 7 (t) — 7 (f; 5). 


V. DOUBLE BELTRAMI FUNCTION 
EXPANSION 


As is shown in the previous section, the pair of vector 
fields (aCi[dMi\, aC e [dM e ]) satisfies the particle rela¬ 
beling symmetry conditions (1151) . This implies that the 
integral path in the configuration space shifted in this 
direction retains the value of the action. Noticing that 
the vector fields pair is obtained by operation of appro¬ 
priate operator on the ion and electron velocities pair as 
follows: 


/ aCildM^ \ _ I Vi \ 

^ aC e [dMe] ) ~ ^ V e ) ’ 


(33) 


This construction procedure is also regarded as an ex¬ 
tension of the general helicity conservation laws found 
by Khesin and Chekanov [20[ to a direct product group 
case. 

Since the exterior derivative of a differential 1-form on 
a three dimensional manifold is given by the curl oper¬ 
ation, substituting m into the mixed helicity (1301) . we 
obtain each part of the mixed helicity as follows: 

( Mi | [dMj]) = —j f (au + a) ■ (aV x u + b) d 3 x, 

a Jx£M 

(Me | [dMe]) = ~2 [ a • bd 3 X. 
a Jx£M 

Note that, in the standard MHD limit a —> 0, an asymp¬ 
totic relation Mi = V) + A ss M e = A = a/a ~ 0(a~ 1 ) 
holds, and thus the hybrid helicity comes close to the 
magnetic helicity as an 0(a~ 2 ) quantity. However, when 
Ci = —C e , the singularity order of H is reduced by a, i.e., 
the leading orders of these helicities cancel each other, 
and the mixed helicity becomes 

H = Ct J d 3 x [2m • b + au ■ (V x m)] , (32) 

which converges to a finite value in the standard MHD 
limit. Thus, the conservation of cross helicity is shown 
to be a special case of the general conservation law for 
the mixed helicity. 


where the integro-differential operator, W, is defined by 


/ Ci(aVx)- 1 + C x aVx -Ci(aVx)- 1 
^ -Ce(aVx)- 1 Ce(aVx)- 1 


(34) 


we recognize that the operation of W on the integral 
path of the action physically implies infinitesimal particle 
relabeling operation (see FigQ]). Thus, we call W helicity- 
based, particle-relabeling operator hereafter. 

It seems reasonable to consider the eigenvalue problem 
of the operator W, because it is expected that the spec¬ 
tral expansion by such eigenfunctions should have some 
“good” properties for the description of the basic for¬ 
mulas and equations. In the following, we will solve the 
eigenvalue problem, demonstrate the mode expansions of 
various quantities and equations, and discuss the relation 
to the uniform background magnetic field effect and the 
standard MHD limit. 

The eigenvalue problem of the operator, 

WV = AV, (35) 

is equivalent to the double Beltrami flow (DBF) problem, 
which is given by the following coupled partial differential 
equations 0 ; 

aV xu + b= -^-m, b = — 7 t-(m - aj). (36) 
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Note that the eigenfunction of the DBF problem is 
constructed using a Beltrami flow, i.e., the eigenfunctions 
of the curl operator, say ?/>( K , a), which satisfies 

V x tp(K,a) = aKxjj(K,a), 


where K, K > 0, and a = ±1 are the mode index, the 
associated eivenvalue, and the helicity of vector field, re¬ 
spectively ||3C|. The Chandrasekhar-Kendall function on 
a cylindrical configuration [31,] and the complex helical 
waves on a periodic box or a Euclidean space [32] are 
known as examples of the Beltrami flows. Expanding 
the variables using xf(K,a), the operator W is reduced 
to a 2x2 matrix as 


/ Ci Ci \ 

a.K +Ciah -ffK 


\ 


Ce_ 

aK 


Ce_ 

aK ) 


for each expansion mode. The eigenvalue of this matrix 
is given by 


\ z rs , a f Ci + C e 

A(A -' T ’ a) = 2l^F- 


CiOtK 


+s 


- CiaK) 2 + 4C?] 2 (37) 


where s = ±1 is the polarity. Note that, the eigenvalues 
for the assigned K and a satisfy 


of the Riemannian metric tensor © for the DBF basis, 
that is to say 


g(K) 


\ 2 

1 aK \ 1 

J(K)~^Cl) + Cf 


(44) 


Substitution of \Fi and \F e into (fT4l) yields the base func¬ 
tions of conjugate momenta space, say \F = \F e ), 

as follows: 



It is easy to check that the base functions SI/ and jF are 
bi-orthogonal each other; i.e., 

(*i(K)\*i(P)) + (*e(K)\* e (P)) 

= 9(1^) 3 K,P ,<?P ^SK ,sp ■ (46) 


Using these base functions, we can rewrite the eigenvalue 
problem (l35l) as follows: 


/ aC.V x <K(K) 
\aC e Vx *e(K) 


A (K) 


*i(K) \ 
*e(K) J ■ 


(47) 


The generalized velocity and momentum are expanded 
using the eigenfunctions of W as 


A (K, a , +)A(AT, a, -) = C l C e , (38) 

and, in the standard MHD limit a —> 0, they become 

A (K, a, +) —> 00 , A (K, a, —) —> 0 for Ci —C e , (39) 
A (K, a, s) —> as\Ci\ for Ci = —C e . (40) 


by 


The eigenfunction of W, say \F = is given 


<F(/\,cr,s) = 


( J 

V A.(K,a,i 


A (K ,<T,s) 


aK 

crC e 


tp(K,a) 


ip(K, cr) 


, ( 41 ) 


V = Y^V{K) ’F(AT), M = J2V( K )fk.( K ), (48) 

K K 

where the expansion coeffcient, V(K), is obtained by the 
inner product, 


V(K) = (V V(K) ). 


(49) 


Since V) = u and aV x M e = qV x (a/a) = —b , the 

expansion coefficients, V ( K , a, s), are determined by the 
following simultaneous equations: 


and we call the set of the eigenfunctions the DBF basis, 
hereafter. If the Beltrami functions if are orthonormal 
each other: 

J ip(K,a K ) ■ ip(P,a P )d 3 x = 5g^6 aKt<7p , (42) 


hereafter, over line and 6 denote complex conjugate and 
Kronecker’s delta, respectively. The corresponding eigen¬ 
functions \F are orthogonal: 


vF(X) *(P))=g(K)6 gp 6 aK , c , P 6 SK , sp , (43) 


hereafter, the tilde denotes the set of mode indices, 
K := ( K,uk,sk ), and g is the inner product of the 
base function \F(A'), i.e., mathematically the component 


u(K , a) = 


1 


aK\ 


A (K, a, +) aC e J 


V(K,a,+) 


1 


~^r) V(K,a,-), 


^A(K,a,~) aC e 
b(K, tr) = ~ (V(K, a, +) + V(K, a, -)) , (50) 


where u(K, a) and b(K, a) are the spectral expansion co¬ 
efficients of the velocity and magnetic fields with respect 
to the basis {fj(K,a)}, respectively. Unfortunately, the 
expansion coefficients V converge to zero or diverge in 
the limit a —> 0, unless Ci = —C e - 

The energy and the mixed helicity are obtained by sub¬ 
stituting (H51) into m and (1501) and using the relations 






















m, ®, and m ; i-e-, 

E = \{ Y V(K)*(K)\Y, V(P)$(P) 

= IY,9(E)\v(K)\ 2 , 

k 

H = aCi ( Y V(K)*i(K) | Y V(P)V x ¥<(P 
+aC e ( Y V(K)*e{K )| Y V(P)V x *e(P 
= Y9(K)A(K) V(K) 


(51) 


(52) 


K 


skew-symmetry of the coefficients of the quadratic terms. 
In the previous study, wherein the generalized Elsasser 
variables (GEV) expansion of the HMHD system was pre¬ 
sented, we conjectured that the conservation of the mod¬ 
ified cross helicity might reflect some symmetry intrinsic 
in the system Q. Since the DBFs become the GEVs for 
Ci = — C e = 1, the modified helicity conservation is now 
recognized as the consequence of a special case of the 
particle relabeling symmetry for ion and electron flows. 


VI. CONSIDERATION OF THE UNIFORM 
BACKGROUND MAGNETIC FIELD 


Applying the DBF expansion to the equation of motion 
yields the following simultaneous equations for the ex¬ 
pansion coefficients, V(K;t): 

^Ufe) = g(K)- 1 YY,(m P \\Q)) 

p Q 

xV{P;t)V(Q-,t), (53) 


For some practical applications, it is important to con¬ 
sider the influence of the uniform background magnetic 
field on the dynamics of plasmas. 

Here, we mathematically consider the influence of a 
background magnetic field, say Bo, which is a harmonic 
function: V x B 0 = 0, V • B 0 = 0. Substitution of 6 + B 0 
into the equations of motion C51) yields 


where the symbol ((if||P||Q)) is given by 

1/1 aK \( 1 aP \ 

Ci\A(K) cr K C e ) V A(P) CrpCe) 


( 1 

aQ N 

1 1 l 

Va(q) 

CQCey 

' C e A(K)A(P)A(Q) _ 


x f ^}{K,a K ) ■ (rp(P,a P ) x t/>(Q,ctq)) d 3 x. 

(54) 

The derivations of (l53l) and (1551) are summarized in Ap¬ 
pendix □ Note that the symbol ((//||P||Q)) is skew- 
symmetric between two arbitrary argument sets, because 
the integrand is given by a scalar triple product of vector¬ 
valued functions whereas the coefficient is symmetric. 
Due to this skew symmetry, we can easily prove the con¬ 
servation laws of the energy and the mixed helicity from 
the expression (1551) as follows: 

f 45 :se(m« a (c) 

K P Q 

xV{K\t) V (P; t) V(Q; t) + c.c. = 0, (55) 

=E E E Pll^ll Q)) mk) a (q) 

K P Q 

xV(K; t) V(P; t) V(Q; t) + c.c. = 0. (56) 

Thus, the DBFs, i.e., the eigenfunctions of the hclicity- 
based particle-relabeling operator are shown to consti¬ 
tute a family of orthogonal function bases that yields a 
remarkably simple spectral representation of the equa¬ 
tion of motion. Especially, the mixed helicity conserva¬ 
tion is naturally built in this representation due to the 


((* 11 ^))=^ 


{dt + Ly i )(a\d Mi] + Bo) — 0, 
(d t + Lv„)(a\dM e ] - B 0 ) = 0. 


(57) 


When the amplitudes of the variables are sufficiently 
small compared to the modulus of Bo, these equations 
are approximated as 

ad t [dAP] = L Bo V t , ad t \dM e ] =-L Bn V e , (58) 

where the identity for two vector fields L^r) = —L^ is 
used (see Eci. (l23l) ~). These linear simultaneous equations 
can be described using the operator W as: 





(59) 


Thus, the linear waves are shown to be expressed by the 
eigenfunctions of the DBF problem (TUI) , with Ci = 1 and 
C e = — 1. The eigenvalue become 


A(A') = — ( aK + s\J {aK ) 2 + 4). 


For typical velocity and magnetic field variables, the lin¬ 
ear simultaneous equations are given by 


d t u = j x Bq - VP, 

(60) 

d t b = V x ((it — aV x b) x .Bo). 

The eigenfunctions of the linearized HMHD system (1501) 
are known as the generalized Elsasser variables (GEV) 
a Physically, the GEVs describe the ion cyclotron or 
whistler waves in plasmas. The phase velocity for an 
assigned K is uj(K) := B 0 k\\ A(Ff) -1 , where fc|j is the 

wavenumber of r/>(A", ctr-) in the direction of Bo- The 
equations of motion for this case are obtained by substi¬ 
tuting V{k-t)e- iul(k)t into V(K;t). 
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It is interesting that the GEV expansion coefficients 
of the basic variables are simply and hierarchically ex¬ 
pressed by V multiplied by the powers of A as follows: 


we obtain the MHD limit of the expansion coefficient, V: 
V(K,ct k ,sk ) = b(K,a K ) - <j k s k u(K,<j k ), (64) 


V e = 


-a\dM e ] = 


u — qV x b 

EkMK)- 1 V{k) il>(k,a K ), 

b 

EkV(K) ip(k,a K ), 


V = u 


(61) 


= EkHK)v(k)ii>(K,<TK), 

a[dMi] = ctV x u + b 

— = EkHK ) 2 v(k) ip(k, <j K ). 

Note that hclicity parameters for the linear wave 
modes have the following significant properties. 

Firstly, the eigenvalues m do not diverge in the stan¬ 
dard MHD limit, a —»• 0: 


. . ~ . Oi OL n 

A(A) « <t k s k + + —<t k s k K 

Z o 


cfrSk- 


This convergence leads to the finiteness of the following 
quantities in that limit: coefficients of the base functions 
©: 


#(A> 


f ( (Jksk + a K ) 

\ {vrSk - fcr K K)tp(K,a K ) 

i/>(K,<jk) 

'i !>{K, (J K ) 

the Riemannian metric (l44l) : 


&kSk 


g(K) = 2 + a 2 K 2 


the coefficient of the RHS of (l5H) : 


/ 1 aK\f 1 aP 

VA(A') &K ) V A(P) cr P 

1 


1 aQ 
A (Q) + VQ 


A(K)A(P)A(Q) 

CFKcrpaQSKSpSQ ( SrK + SpP + SqQ 
' SpP + SqQ 




-a a K crpcrQKPQ 
| sqQ + srK 


s k K 

sqQ + srK skK + SpP 


+ 2 ) + o(cc 2 ) 


SpP SqQ 

CTKVpVQSKSpSQ I SrK + SpP + SqQ ) . 


(62) 


Since the simultaneous equations (l50l) converge to 

u(K,ct k ) = a K s K V(K\crK,-) ~ a K s K V(K, a K , +), 
b(k,a K ) = V(k,cr K ,+) + V(k 1 a Kl -), (63) 


and thus, the equation of motion (15^1) has the standard 
MHD limit. The coefficients, V, are associated with the 
conventional Elsasser variables by the formula 

z + = u + b 

= V(K,+,-)ip(K,+) + V(K,-,+)'i/)(K,-), , 
z- = u — b 

= —V(K, +, +)ip(K, +) - V(K, -, -W(K, -). 

The base function in momentum space >&, on the other 
hand, diverges on the order of a -1 , at which the diverging 
a/a term is reflected. 

Secondly, the singularity order of the mixed helicity 
(15011 reduces by a and the constant H become the mod¬ 
ified cross hclicity, which converges to the cross helicity: 

Hc= J d 3 x [2 u ■ b + au ■ (V x it)] 

d 3 x u ■ b. (66) 

In our previous study, it was shown that the conservation 
of the modified cross helicity is naturally derived from the 
GEV representation of the HMHD dynamics. 



VII. DISCUSSION 

In the present study, we considered the hclicity con¬ 
servation laws of HMHD system from Lagrangian me¬ 
chanical, invariant action theory viewpoint. The hybrid 
and magnetic helicity conservation laws were derived as 
consequences of the particle relabeling symmetry of the 
ion and electron flows, respectively. To prove the conser¬ 
vation laws, it is convenient to use the pair of ion and 
electron velocity fields © as basic variables, while that 
of fluid velocity and current fields had been used in our 
previous study [U. Mathematically, this variables change 
was carried out by changing the configuration space of 
HMHD system from semidirect product group to direct 
product one. 

Furthermore, associated integral path variation of the 
invariant action was shown to be expressed by the op¬ 
eration of the helicity-based, particle-relabeling opera¬ 
tor (1341) on the reference path, which maps the general¬ 
ized velocities (V), V e ) to the action-preserving, particle¬ 
relabeling fields (aCj\ dMj \ , aC e \dM e ]). 

The eigenfunctions of the relabeling operator are 
DBFs, which are well-known, force-free solutions of the 
HMHD system m. and found to provide a family of 
orthogonal function bases that yields the spectral rep¬ 
resentation of the equation of motion with a remark¬ 
ably simple form. Thus, the GEV based formulation we 
had discussed in [|J is now understood as an example of 
more wider class of orthogonal function expansion of the 
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HMHD equations, since the GEVs are special case of the 
DBFs (Ci = -C e = 1). 

The implication of this eigenvalue problem may be 
well-understood by considering the correspondence be¬ 
tween the Lagrangian and Hamiltonian mechanics. 

It is well-known that the Lie algebraic structure nat¬ 
urally induces a so-called Lie-Poisson structure on the 
dual space of the Lie algebra by defining the Poisson 
bracket by {A,B}(n) = (m|[§^, f§]), where p € 0 * 
is an element of the generalized momentum space, and 
A, B £ T(g*) are the functionals of the generalized mo¬ 
menta [jjj. In the incompressible HMHD case, the Pois¬ 
son bracket based on © and 0 becomes 


{A,B}(M) 


M, 


V x 


SA SB » \ 

Zm~ x JmJ) 


+ M e 


V x 


SA 

SM e 



(67) 


When B is the Hamiltonian obtained by the Legendre 
transformation of the Lagrangian m, which results in 
B(M) = L(V), the functional derivatives of B are given 
by JM 7 = Vi > Ti (fl = Ve ' Integration by parts of 
yields 


{ A,B}(M ) = J d 3 x 


+ 


6A (Vi x (V x A*)) 

■ (V e x (V x M e )) 


SM t 
SA 


5M e 


.( 68 ) 


By setting the derivatives = £, = V, the cubic 

terms of the first variation (fl6l) are reproduced, and thus, 
the action-preserving variation is shown to correspond to 
the functional derivative of a certain Casimir function. 

In the Hamiltonian mechanical approach to the sta¬ 
bility problem of the equilibrium solutions, the DBF are 
known to constitute the dynamically accessible variations 
that a priori satisfy the conservation laws for energy 
and Casimirs Q. In the incompressible HMHD case, 
the Casimirs are given by the magnetic and hybrid he- 
licities. In the Lagrangian mechanical approach, on the 
other hand, they are obtained from the invariant action. 
Thus, the eigenvalue problem for the invariant action is 
naturally described as the DBF problem. 

Since the DBFs for assigned C* and C e were orthogo¬ 
nal each other, by using them as base functions we could 
obtain a general form of the “normal mode” expansion 
of the Riemannian metric, the structure constants of the 
Lie algebra, and the equation of motion. The combina¬ 
tions of Ci and C e are arbitrary, i.e., the DBF basis has 
two degrees of freedom. By changing the values of Ci and 
C e , we obtained a family of “canonical” transformations 
between the spectral representations of the equation of 
motion. The spectral representations of the equation of 
motion formally have a common mathematical expres¬ 
sion given by m . which is known as the Euler-Poincare 
equation (Chapter 13 of Ref, [hi 1. or as the geodesic equa¬ 


tion fioj : 


±JL = r b 

d t dC l da ? d£ b ’ 

where the variables and coefficients , C\ a 

respectively correspond to V(K), M{K), and 

((K\\P\\Q))A(K)/g(K) in the present study. As is 
expected, conservation laws for the energy and the 
mixed helicity are easily proved from the symmetric 
properties of the obtained structure constant. 

In the standard MHD limit a —»• 0, the eigenvalues, 
and thus, the related quantities such as the expansion 
coefficients, the Riemannian metric, and the structure 
constants, diverge or shrink to zero unless Ci = —C e . 
Hence, the DBF basis seems unsuitable for comparative 
analysis between HMHD and MHD in most cases. How¬ 
ever, it is very interesting that consideration of the effect 
of a uniform background magnetic field yields such a lin¬ 
ear wave equation that uses the same DBF operator as 
Ci = — C e = 1. It is well-known that the linear wave 
modes in the incompressible HMHD system are the ion 
cyclotron and whistler waves and are elegantly described 
by the GEV [2lJ. That is, the GEV is such that the 
DBF has non-diverging properties in the limit a —>■ 0. 
Thus, among the wide variety of the DBF expansions, 
the GEV expansion is the most suitable for comparing 
the dynamics of the HMHD system to its MHD limit, 
avoiding singularity problems. 

Since the DBFs are constructed from the eigenfunc¬ 
tions of the curl operator, it is easy to include Laplacian- 
type dissipation into the spectral representation of the 
equation of motion; 

{g{K)^ + D (K,v, rj))v{K-,t) 

= EE {(K\\P\\Q))A(Q)V(P-,t)V{Q-,t), (69) 

p Q 


where D 1 v, and iq are the dissipation term coefficient 
given by 


D(K,v,ti) := K 2 



( t k K \ 2 77 

Ce) Cl\ ’ 

(70) 


the kinematic viscosity, and the resistivity, respectively. 
This feature allows us to apply the DBF expansion to 
such analyses as the closure problem [ 2 lJ or the direct 
numerical simulation ( 22 ]. 
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Appendix A: Derivation of the Lin constraints 


Appendix B: local expression of particle-relabeling 
symmetry 

By the term “particle-relabeling symmetry,” we rec¬ 
ognize the invariance of the flow against the change of 
Lagrangian coordinates. The freedom of choice of the 
action-preserving transformation exists only at the “ini¬ 
tial time” and the transformation along the integral path 
of the action is determined by this initial condition. This 
symmetry is qualitatively different from the symmetry 
considered, for example, in gauge field theory, wherein, 
in principle, group transformation is applicable at any 
point in the relevant space and time [331 ]. 

Thus, the evolution of transformation should be con¬ 
sidered. Let £ and e be a displacement-generating vector 
field and a small parameter, respectively. For an assigned 
flow and a sufficiently small displacement, the displace¬ 
ment field must satisfy the particle tracing relation: 

X(a + e£(a;0);t) = X(a;t) + e£(J?(o;i);i), (Bl) 


We briefly review the derivation of the formula for the 
variation of the tangent vector to an integral path. For 
this derivation process, we use only the exponential map 
and the Baker-Campbell-Hausdorff formula. Since no 
material specific to any particular Lie algebra is used 
here, the result is applicable to all the Lie groups. 

Let 7 (t; 6) (t £ [0,1]) be a path on G with a variation 
parameter 8 £ I £ 1. The path C—xD is approximated 

by 


7 (t + r-,6) ss exp [r(V(f) +<5V'(t))] o 7 (£; 5), 

where V (t) is the tangent vector to the reference path 
(<5 = 0), and SV(t) is the small deviation. The path 
C—xA—xB—xD is, on the other hand, also approximated 
by 


where the X is the PTM for the assigned flow. At the 
order O(e), each component of £ satisfies 

= e(X(a;t);t). (B2) 

(o;t) 

Differentiating with respect to t and evaluating at t = 0, 
we obtain 


£ fc (a;0) 


dX l 

dx k 


e(a; 0 ) 


d 2 X i 

dG 

dX k 

dx k dt 

(o;0) dxk 

(X(a- 0);0) 


dt (X (a; 0 ); 0 )" 


(B3) 


The relation (JoJ) leads to the following PDE for the vector 
fields in the Eulerian specification: 


7 (t + r; 8) ~ exp (d£(t + r)) o exp (rV(t)) 
oexp(-6£(t)) 07 (t;S) 


dC_ 

~dt 


vk d^_ i kd_V^\ 

dx k dx k ) 


= 0, (B4) 


(X(a;0);0) 


(see Figure O. Expanding £(t + r) = £(t) + Tdt£(t)+o(T) 
and using the Baker-Campbell-Hausdorff formula at the 
lowest two orders, we obtain 


exp (8£(t + t)) o exp (rV (t)) o exp ( — 8£(t)) 


= exp 


TV(t)+T6(d t £(t) + [£(t), V(t)]) +o(5) +o(r) . 

(Al) 


which is the evolution equation for the frozen-in line ele¬ 
ment. Since the Lie bracket of the vector fields is given by 
(155b . the obtained evolution equation agrees with (IA2I) 
for V(t) = 0. For divergence-free fields in a three- 
dimensional space, the equation is rewritten using vector 
analysis notation as 

^ + Vx(CxF) = 0. (B5) 


Since the two approximated paths from C to D agree 
with each other in the limit S —> 0 and r —x 0 , we obtain 
the Lin constraints 


Appendix C: Structure constants for the DBF basis 


V(t) = d t £(t) + [£(t), V(t)] (A2) The symbol ((A'||P||Q)), which is related to the struc¬ 

ture constant of the Lie algebra, is defined by using a 
combination of the Riemannian metric and the Poisson 


at the order O(St). 
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bracket for the HMHD system as follows: 

((K\\P\\Q)) A(K) := {^(if)|[#(P),#(Q)]; 

= (2i(A-)|vx ('S'i(P) x 

-(^e(it:)|vx (* e (P) X * e 

d 3 4 5 6 2 7 8 9 10 f [cr 1 ^^) • (^(P) X 'Jl(Q)) 

+C~ 1 * e (K) ■ (^e(P) x ¥ e (Q))l, (Cl) 


+ ( 
A(PT) 


We obtain the Euler-Lagrange equation (l53l) if the fixed 
path end conditions are imposed. 


where the second line is derived by integrating by parts 
with respect to x and by using the relation (|47| . Sub¬ 
stitution of m into the second line yields (f54|) . i.e., the 
explicit expression of the symbol ((AT||P||(3)). Using this 
symbol, the first variation of action reads as 


dt( Y V(Q ) dttik) + [£{K),Y V(P)] 


dt(g{K)V{K)d£{K) 


+ Y HQ){(Q\\K\\P))Vm(K)V(P )), 

p,Q 


= (g(K)V(F)&K)) t=1 - (g(K)V(F)aK)) t=o 

+ j dt£(K)(- g(K)d t 9(F) 

+ Y MQ){{K\\P\\Q))V(Q)V(Pj) . (C2) 
p,Q 
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